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Abstract. The nature of the Hyades stream, or Hyades moving group, is a long-standing question of Galactic Astronomy. 
While it has become widely recognized that the Hercules stream, an unbound group of stars lagging behind galactic rotation 
and moving outward in the galactic disk, is associated with the outer Lindblad resonance of the rotating galactic bar, there 
is still some debate about the nature of the more prominent low-velocity stream sharing the kinematics of the Hyades open 
cluster. Is this stream caused by additional non-axisymmetric perturbations of the galactic potential, such as transient or quasi- 
stationary spiral waves, or by the on-going evaporation of the Hyades cluster? A simple observational test has been designed to 
answer that question, i.e. to determine whether the Hyades stream is primarily composed of coeval stars originating from the 
Hyades cluster, or of field stars. Using the Geneva-Copenhagen survey of F and G dwarfs, we compare the mass distribution 
and metallicity of the stream to those of field disk stars. If the Hyades stream is composed of stars trapped at resonance, its 
mass distribution should obey the present-day mass function (PDMF) of the disk, and its metallicity should reflect its origin 
in the inner regions of the Galaxy. On the other hand, if it is an evaporated cluster, we expect a different mass distribution, 
depending on the inital mass function (IMF) of the cluster, and on the proportion of evaporated stars as a function of mass. 
We find that extreme conditions have to be adopted for the selective evaporation and IMF of the cluster to make the observed 
mass ditribution of the stream only roughly consistent (at a one-sigma level) with the coeval evaporated cluster scenario. The 
observed mass distribution is in much better agreement with the PDMF of the field. We also note that the peculiar metallicity of 
the stream is inconsistent with that of a field population from the solar neighbourhood trapped in the primordial cluster during 
its formation process and subsequently evaporated. These observations thus favour a resonant origin for the Hyades stream, as 
suggested in Famaey et al. (2005). 

Key words. Galaxy: kinematics and dynamics - evolution - clusters and associations - disk - solar neighbourhood - stars: 
kinematics 



1. Introduction 

It has been known for a very long time that a spatially un- 
bound group of stars in the solar neighbourhood is sharing 
the same kinematics as the Hyades open cluster (Hertzsprung 
1909, Stromberg 1922, Eggen 1958, Perryman et al. 1998). 
Assuming that it was a vestige of an initially more massive 
cluster which partly evaporated with time, Eggen christened 
this kinematically cold group the Hyades supercluster. More 
generally, it is called the Hyades stream, or Hyades moving 
group. 

During the last fifteen years, Eggen's hypothesis that kine- 
matic groups of this type are in fact cluster remnants has been 
largely debated, because they may also be generated by a num- 
ber of global dynamical mechanisms. A rotating bar at the cen- 
tre of the Milky Way may e.g. cause the velocity distribution 
in the vicinity of the outer Lindblad resonance to become bi- 
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modal, due to the coexistence of orbits elongated along and 
perpendicular to the bar's major axis. Today, this mechanism 
is thought to actually account for the Hercules stream (Dehnen 
2000), a group of stars lagging behind the galactic rotation and 
moving outward in the disk. Moreover, the perturbation that the 
triaxial bar induces on a flat axisymmetric disk does produce 
some chaos: Fux (2001) showed that when the bar is taken into 
account, the chaotic regions, decoupled from the regular re- 
gions, are more heavily crowded in the region of the Hercules 
stream in velocity space. On the other hand, other streams, in- 
cluding the Hyades, could still be linked with the spirality of 
the Galaxy, because any perturbation of the axisymmetric po- 
tential is likely to buffet the stars. 

Eggen's scenario and the dynamical perturbations are, of 
course, not mutually incompatible phenomena. Clusters are 
known to evaporate over time, and there must be some inter- 
mediate state when a group of stars no longer spatially identi- 
fiable still share similar velocities. Indeed, disk stars (most of 
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which move on quasi-circular epicyclic orbits) which formed at 
the same place and time, and which stayed in the same region 
of the Galaxy after a few galactic rotations must necessarily 
have the same period of revolution around the Galactic cen- 
ter, and thus the same guiding-center which in turn implies the 
same tangential velocity (Woolley 1961). On the other hand, 
our Galaxy is known to have large spiral perturbations that are 
likely to have a kinematic effect on the velocity of stars. In the 
past, the importance of stirring by spiral structure has been un- 
derestimated because it was thought that spirals heated the disk 
strongly, and because the amount of heating is observationally 
constrained. However, Sellwood & Binney (2002) showed that 
the dominant effect of spirals is to stirr without heating. Thus, 
the Hyades stream may well be an outward-moving stream of 
stars on horseshoe orbits that cross the corotation of the spi- 
ral pattern. The kinematics of K and M giant stars in the solar 
neighbourhood seem to comfort this scenario, and suggest that 
the Hyades stream is composed of metal-rich stars with a wide 
range of ages (Famaey et al. 2005, see also Chereul & Grenon 
2001). However, metallicities were available for too few giants, 
and individual age estimates were not precise enough to reach a 
decisive conclusion. This motivates the present analysis, where 
we examine whether properties of F and G dwarfs from the 
Geneva-Copenhagen survey (Nordstrom et al. 2004) are in bet- 
ter agreement with this dynamical scenario, with Eggen's sce- 
nario, or with a mix of both. 



2. The sample 

The Geneva-Copenhagen catalogue of F and G dwarfs 
(Nordstrom et al. 2004) is the result of a decade-long cam- 
paign of Stromgren photometric and Coravel spectroscopic 
measurements for more than 16000 late-type dwarfs in the solar 
neighbourhood. It gives positions, parallaxes, proper motions, 
masses and metallicities for a sample of dwarfs complete to 
40-70 pc depending on spectral type. Metallicities are deter- 
mined from Stromgren photometry according to a tailored cal- 
ibration described in Nordstrom et al. (2004; their Sect. 4.3). 
Masses are determined by comparing the observed position in 
the Hertzsprung-Russell diagram (effective temperatures and 
absolute magnitudes were derived from Stromgren photom- 
etry and Hipparcos parallaxes) with stellar-evolution tracks 
from Girardi et al. (2000). Masses are displayed as a func- 
tion of the b — y color index in Fig. [2 which shows that the 
Geneva-Copenhagen survey is only complete in the mass range 
[0.8, 1.5] M . 

The survey contains 122 stars flagged as definite or possi- 
ble members of the Hyades cluster by Perryman et al. (1998) 
and de Bruijne et al. (2001). To avoid contamination from the 
Hyades cluster in our subsequent analysis of the Hyades stream 
we exclude those 122 stars from the sample. After also exclud- 
ing the binaries to avoid the distorting effect of binarity on the 
photometry, the total number of (single) stars in the present 
sample is 8084. 
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Fig. 2. Isocontours for the whole survey in the t/ V-plane: the 
contours correspond respectively to 0.5, 0.8, 1.2, 1.5, 1.9, 2.6, 
3.1, 3.5, 3.8, 4.2, 4.7, 5. stars/(km/s) 2 . The Hyades stream, at 
U ^ -37 km/s and V - -17 km/s, is the most prominent fea- 
ture. 

3. The velocity field 

Fig.|2]presents the distribution of all stars from the survey in the 
t/ V-plane (U is the velocity towards the galactic center, V the 
velocity in the direction of Galactic rotation, both with respect 
to the Sun). For computing U and V, distances are adopted 
from the Geneva-Copenhagen Catalogue. These distances are 
derived from Hipparcos parallaxes if their relative errors are 
better than 13%, and from photometric distances (which are 
uncertain by 13%) otherwise. 

The small-scale structure in the t/V-plane appears clearly 
on Fig. 2: the main overdensities are the well-known Hercules 
stream (U - -30 km/s, V ^ -50 km/s), the Sirius-UMa 
stream (U - 10 km/s, V - -5 km/s), the Castor group (U ^ 
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-10 km/s, V — -10 km/s), the Pleiades stream ( U - -15 km/s, 
V = -25 km/s), and, the most prominent, the Hyades stream 
(U -37 km/s, V - -17 km/s). Those well-known overdensi- 
ties (see Dehnen 1998, Chereul et al. 1999, Montes et al. 2001, 
Famaey et al. 2005) are either due to the evaporation of star 
clusters, or to non-axisymmetric perturbations by spiral arms 
and the galactic bar (e.g. the Hercules stream). Their promi- 
nence close to the center of the i/V-plane is not compatible 
with the hypothesis that they are the result of merger events. 
Substantial amounts of such satellite debris have been identi- 
fied in the Geneva-Copenhagen survey (Helmi et al. 2006), but 
at higher velocities and with lower metallicities. 

4. The Hyades stream 

We shall now focus our attention on the prominent Hyades 
stream, which was suggested by Famaey et al. (2005) to be 
a dynamical stream coming from the inner Galaxy because of 
a spiral perturbation, although firm proof was lacking because 
metallicities were available for only very few stars and indi- 
vidual age estimates were not precise. The main question to be 
answered here is thus whether this kinematic group is really a 
dynamical stream, or rather a coeval evaporated cluster. 

An evaporated cluster would have a cluster-like homogene- 
ity of age and metallicity, while a dynamical stream would have 
the age and metallicity composition of a random disk popula- 
tion (its metallicity could however reflect that of the galacto- 
centric radius wherefrom it originates as a result of the dynam- 
ical perturbation). 

However, even with the high-accuracy data of the Geneva- 
Copenhagen survey, the precision on the metallicity and age is 
not high enough for one of the scenarios to stand out immedi- 
ately. Indeed, the ages determined in Nordstrom et al. (2004) 
from the position in the HR diagram cannot be reliably used 
(Pont & Eyer 2004): the uncertainties and systematic biases on 
the age determination are strongly correlated with the mass, 
and a group of stars with different masses and homogeneous 
age cannot be detected as such from the computed ages (see 
Fig. 25 of Nordstrom et al. 2004). 

Nevertheless, the two scenarios can be distinguished by a 
more subtle method than the simple analysis of age and metal- 
licity distributions. If the Hyades stream is caused by a non- 
axisymmetric dynamical perturbation of the potential, we ex- 
pect a mass distribution similar to the field stars, obeying the 
present day mass function (PDMF). On the other hand, if the 
Hyades stream is an evaporated cluster, composed of coeval 
stars (600 Myrs old), we expect a different mass distribution, 
depending on the inital mass function (IMF) of the cluster, and 
on the proportion of escaping stars as a function of mass (see 
e.g. Terlevich 1987). 

4.1. Predictions in the purely coeval case 

The mass distribution thus offers us a good test of coevality for 
the Hyades stream. Let us assume that the stream is coeval, en- 
tirely generated by the on-going evaporation of the primordial 
Hyades cluster. We consider two groups of stars, chosen to be 



as far apart as possible without suffering from selection biases 
and incompleteness (as illustrated by Fig.^: 

- Group 1: stars with 1.3 M < m < 1.5 M . The upper limit 
is below the Hyades cluster turnoff mass (2.2 M ), so that 
no giant star from the assumed coeval stream is present in 
the group. To avoid contamination by the small number 
of field giants present in the Geneva-Copenhagen survey, 
we further restrict this group to b — y < 0.42 (this sup- 
plementary condition does not alter the selection for main- 
sequence stars, see Fig.l). 

- Group 2: stars with 0.8 M < m < 1 M . 

Then, in this coeval scenario, the proportion of stars belong- 
ing to the Hyades stream with respect to field stars of the disk, 
Miya/^Vfieid, should decrease from the value p\ in Group 1 to 
the value p2 in Group 2, according to the variation of the ra- 
tio of the mass function of the young evaporated Hyades clus- 
ter to the PDMF of field stars. The estimation of the ratio 
Pi/p2 under this coeval hypothesis nevertheless depends on the 
assumptions made on the efficiency of selective evaporation, 
namely on the proportion of evaporated stars from the clus- 
ter in the mass range corresponding to Groups 1 and 2. It is 
well-established that low-mass stars are preferentially depleted 
from star clusters from the effect of mass segregation (see e.g. 
Baumgardt & Makino 2003), leading to a flattening of the ob- 
served cluster mass function, or even turning an initially in- 
creasing mass function into one which is decreasing towards 
low-mass stars (see e.g. Reid & Hawley 1999, Dobbie et al. 
2002). However, it is not at all obvious that, for a 600 Myrs 
open cluster, a significant preferential depletion should be ex- 
pected for stars in the range [0.8, 1] M Q as compared to stars in 
the range [1.3, 1.5] M . 

4.1 .1 . Without selective evaporation 

Let us first assume that the relative depletion of stars from the 
primordial Hyades cluster is the same for the two mass groups 
under consideration. The values p\ and pi are then entirely de- 
termined by the ratio of the Hyades cluster IMF to the field 
PDMF, integrated over the appropriate mass intervals. 

Moreover, since by construction Group 1 is restricted to 
dwarf stars, it is also necessary to correct the PDMF for the ab- 
sence of field giants when estimating p\\ the ratio of the main 
sequence lifetime over total lifetime is 80% for a 1 .5M star 
with Z = 0.008 (Lejeune & Schaerer 2001, assuming that the 
evolution posterior to the Helium flash represents 10% of the 
MS lifetime), implying a correction of the PDMF by a factor 
0.8 in Eq.(l) below. Such a correction is unnecessary for Group 
2, dominated by low-mass stars that merely have the opportu- 
nity to leave the main-sequence when they are as old as the 
Galaxy itself. 

Adopting for the PDMF dN /dm oc m -45 (Kroupa et al. 
1993), and for the present-day star formation IMF dN/dm oc 

m -23±0.7 foj . L3 Mq < m < l 5 M& and dN f dm a- OT -2-7±0.3 f or 

0.8 M <m<\ M (Kroupa 2001), we have: 
f L5M0 IMFdm f 1M ° PDMF dm 

Pl_ _ JhWy liVlrU ' W J().8M rL,Wlr Um = 18+1 . 03 

P* 0.8 PDMF dm C IMF dm ' ^ 
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4.1 .2. With selective evaporation 

N-body simulations of the dynamical evolution of open clusters 
with an initial number of 1000 stars (Terlevich 1987) have long 
ago demonstrated that evaporation was more effective for low- 
mass stars because of mass segregation. This trend is extremely 
clear in numerical simulations if one compares the relative de- 
pletion of stars lighter and heavier than m ~ 0.5 M Q . However, 
if one compares the relative number of escapers from our two 
mass groups after 75% of the stars have left the cluster, the ef- 
fect is not clear at all (see Fig. 2 of Terlevich 1987). The numer- 
ical noise is high, and the depletion is favoured in one group or 
the other depending on the initial conditions. The simulations 
seem to indicate that the relative number of escapers in Group 
2 could be higher than in Group 1 at the very most by a factor 
k — 1 .4. While such an efficient selective evaporation would 
steepen the observed mass function of stars evaporated from the 
cluster to bring it closer to the PDMF (and would decrease the 
predicted ratio p\ / p 2 by a factor k), it would actually &\soflat- 
ten the observed mass function of the cluster itself. From Fig. 
8 of Reid & Hawley (1999), it appears that (for the mass range 
under consideration here), the Hyades cluster MF has a power- 
law exponent ~ -2.3. A factor k for the relative number of esca- 
pers in the range -0.1 < log(m) < as compared to the relative 
number of escapers in the range 0. 1 < log(m) < 0.2 would then 
lead to an IMF with a power-law exponent of -2.3 - 5 log(/c). 
An extreme scenario, combining a selective evaporation with 
a maximal k — 1.4 and the steepest IMF oc m~ 10 (allowed 
by Rroupa 2001 error bars) is thus not formally excluded. In 
that extreme case, the lower bound in Eq.(l) (corresponding 
to an IMF power-law exponent of -3.0) should be divided by 
k — 1.4, thus yielding p\lp2 = 1.76. 

4.2. Observational test of coevality 

Let us now determine the observed value of pi and p2 from the 
Geneva-Copenhagen survey. In order to estimate Nhya/Nmd in 
the two mass groups, we need to separate the Hyades stream 
from the field of the galactic disk. We therefore consider the 
strip -21 < V(km/s) < -12 in the f/V-plane, corresponding to 
the range in V covered by the Hyades stream (see Fig. 2). We 
only consider stars with [Fe/H] > -0.5 to avoid contamination 
by the halo (the Hyades overdensity is not seen in the {/V-plane 
for stars with [Fe/H] < -0.5). 

Figs. 3 and 4 display the histogram along U in the V- 
strip for Groups 1 and 2, respectively. The relative Poisson 
noise is larger for Group 2 due to the smaller number of stars. 
Neglecting the Castor moving group appearing as a sharp peak 
in the f/-bin [-15,-5] in Group 2, these histograms may be 
represented by a superposition of two normalized gaussians 
N({U), erf/), the first one corresponding to the Hyades and the 
second one to the field of the galactic disk: 

#i(*i) = #i [(1 -*i)xJV(-37,12) + ifci xM-2,22.5)] (2) 

for Group 1 , and 

N 2 (k 2 ) = N 2 [ (1 - h) x M-37, 12) + k 2 x N(-2, 32) ] (3) 

for Group 2, where N\,N 2 are the total number of stars of the 
considered group in the V-strip, multiplied by the width of the 
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Fig. 3. Histogram of the U velocities for the 415 stars in Group 
1, with -21 < V(km/s) < -12. The dashed line corresponds to 
the contribution of the Hyades stream. 
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Fig. 4. Same as Fig.[5]for the 323 stars in Group 2. The sharp 
excess in the t/-bin [-15,-5] is due to the presence of the 
Castor moving group, and to a high relative Poisson noise. 



U -bins, and k\ , k 2 the fraction belonging to the field. The mean 
value (£/) = -2 km/s for field stars does correspond to the 
value found for giant stars in Famaey et al. (2005) when the 
streams were not taken into account. The velocity dispersion 
of the field increases with decreasing mass (cr^ = 22.5 km/s 
in Group 1 and o~u = 32 km/s in Group 2), as expected from 
the age-velocity dispersion relation (see e.g. Dehnen & Binney 
1998). 
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The parameters k\ and ki have been derived numerically 
(separately for Group 1 and 2: j = 1,2) by minimizing: 



[Ni-Niikj) ] 2 



Ni(kj) 



(4) 



where the sum extends over the M bins in U, Nj is the observed 
number of stars in bin i, and the denominator N, corresponds 
to the expected number. In Group 2, the bin containing Castor 
has not been included in the sum. Pearson's x 1 goodness-of-fit 
criterion used here should only be applied to bins containing at 
least 5 stars. Therefore, the fitting process has been restricted 
to the range -65 < £/(km/s) < 55. 

The uncertainties cr^ , cr^ on the parameters k\ , ki (listed in 
Tabled have been estimated from the relation 



cr, 



= s * N . I^i) 2 



(5) 



where the derivatives dkj/dNj are derived from the condition 
dx 2 /dkj = 0, yielding 



2 



Ni 



M-37,12) 



M-2,22.5) - 7V(-37,12) 



V 



(6) 



under the assumption that the number of stars in the bins are 
large enough to ensure that \Nj - Ni\/Nj « 1. 

In order to derive p\ and p2, we need to estimate the propor- 
tion of field stars in the considered strip -21 < V(km/s) < -12 
with respect to the whole £/V-plane. To avoid contamination 
by other streams present in our sample, we simply approxi- 
mate the field-star V-distribution by N(-10, 12) for Group 1 
and Af(-15, 18) for Group 2, according to the age-velocity dis- 
persion relation and the asymmetric drift (Dehnen & Binney 
1998). We thus obtain 



f 2 

J-21 



pi = (1 -h)lh X | /V(-10, 12)dV = 0.175 ±0.024, (7) 
and 

p 2 = (l-k 2 yk 2 x f /V(-15,18)dV = 0.136 ±0.025, (8) 

J-21 

where the errors have been propagated according to cr p = 
cr k \dp/dk\ = a-k/k 1 £^NdV. We thus have: 



E± _ 1 9Q+0.50 
„ - i " ZV -0.35- 
Pi 



(9) 



In order to evaluate the sensitivity of these results to the choice 
of the V strip ([—21, -12] so far), Table^lists k\,k%,p\,pi for 
slightly different choices of this strip. All these choices actually 
yield very similar results for the p\jpi ratio, clearly incompati- 
ble with the standard coeval scenario without selective evapora- 
tion, predicting 3.18*i'S [see Eq.(l)]. Even the most favourable 
coeval scenario, with a steep IMF of power-law exponent -3.0, 
and with a significant preferential evaporation of 1 M stars 
over 1.5 M stars, yields p\/pi = 1.76, a ratio higher than the 
observed value by one sigma. In this scenario, the mass of the 
stream (around 800 M in the Geneva-Copenhagen catalogue 



alone, i.e. a lower bound since it contains only FG dwarfs) im- 
plies a very massive primordial Hyades cluster. Chumak et al. 
(2005) showed that, after 600 Myrs, 800 M could be lost by 
the cluster if the initial mass was 1400 M and the initial virial 
radius 7.5 pc: this evaporated mass would be distributed in a 
volume of radius 300 pc around the Sun, but since 800 M 
is observationally a lower bound, even more extreme initial 
conditions would probably be needed to account for the full 
Hyades stream. We thus conclude that the most likely sce- 
nario is that the Hyades overdensity in the {/V-plane is not 
uniquely composed of coeval stars evaporated from the primor- 
dial Hyades cluster. 

Note that the zeroth order expectation p\ - p2 for a stream 
exclusively composed of field stars is well within the error bars 
of Eq.(|9). However, the best fit value yields p\ > p 2 - There 
might be several reasons for this: 

- First of all, since the ratio of Eq.(|9j was observationally es- 
timated under the coevality hypothesis, we did not take into 
account possible variations of the intrinsic velocity disper- 
sion of the Hyades stream itself. However, under the dy- 
namical perturbation hypothesis, the Hyades stream itself 
could be affected by the age-velocity dispersion relation. 
For instance, if we take for the Hyades <Ju = 15 km/s 
instead of 12 km/s in the older Group 2 (Eq|3), we find 
P\IP2 — 1-15^0 32> a value closer to one than Eq.J5J. 

- pi > Pi might be caused by the PDMF of the Hyades 
stream being slightly different from the PDMF of the so- 
lar neighbourhood. As we shall see in the next section, the 
Hyades stream is over metal-rich at all masses, pointing 
towards an origin in the inner Galaxy. Since the lifetime 
of stars depends on the metallicity, the PDMF of the inner 
Galaxy should be slightly different from the PDMF of the 
solar neighbourhood, even with identical IMF and stellar 
formation rates (SFR). Moreover, the SFR could also vary 
as a function of the galactocentric radius. 

- P\ > Pi might be due to some of the stars in the stream still 
belonging to the evaporated cluster since we know that the 
Hyades cluster is kinematically associated with the Hyades 
stream, and we also know that this cluster does evaporate 
with time and must create a velocity clump of spatially ex- 
tended stars. Assuming the same PDMF for the Hyades 
stream and the solar neighbourhood (zeroth order approx- 
imation), and the best fit value for the ratio p\/pi, we can 
actually estimate what proportion of stars in each group is 
associated with the coeval Hyades cluster, and what propor- 
tion is associated with the dynamical stream. Eq. (1) gives 
the ratio of the evaporated cluster proportion in Group 1 
with respect to this proportion in Group 2, while the ratio 
is assumed to be 1 for the dynamical stream. Eqs. and 
{8} then give the sum of these proportions in each group. 
From there we find that the evaporated cluster represents 
about 40% of the total stream in Group 1, but only 15% 
in Group 2. These numbers are given here as an indication 
since these proportions should be affected by the variation 
of the stream PDMF, and could of course be zero within the 
error bars of Eq.(|9j. 
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Table 1. Sensitivity of the parameters characterizing the relative fraction of Hyades and field stars to the choice of the V strip 
enclosing the Hyades (see Eqs. |2] [5] |4] and [8] for the definition of the various parameters listed). For Group 1, Pearson's 
goodness-of-fit^ 2 (A:i) behaves roughly like a% 2 distribution with 10 degrees of freedom (12 bins - 2 constraints). For Group 2, 
that number reduces to 9 (because of the exclusion of Castor's bin). The probability that a random variable exceeds the observed 
X 2 value is in all cases large enough for the two-gaussion model to be considered as a satisfactory fit to the data [To fix the ideas, 
ProbCf^ > 17.4) = 6.6% and Prob(^ > 13.7) = 13.4%]. 
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4.3. Metallicity of the stream 

We have thus shown that the vast majority of stars in the 
Hyades stream are most likely field-like disk stars obeying the 
PDMF of the disk. This does not necessarily imply that the 
stream originates from a non-axisymmetric perturbation of the 
Galactic potential. One could instead imagine that the primor- 
dial Hyades open cluster was an extremely massive object able 
to trap some older galactic field stars during its formation pro- 
cess (see e.g. Fellhauer et al. 2006). These stars could have 
formed a dynamically hotter sub-system in the outer part of the 
object, and could have been preferentially evaporated (whilst 
obeying the PDMF). However, in that case, one would expect 
the stream to exhibit a slightly sub-solar metallicity, charac- 
teristic of the solar galactocentric radius, close to which the 
Hyades cluster is supposed to have formed if it has not been 
dynamically perturbed by spiral arms. 

The mean metallicity of the Geneva-Copenhagen survey is 
[Fe/H] = -0.16. Excluding halo stars with [Fe/H] < -0.5, 
and excluding all the stars from the Hyades box (-21 km/s < 
V < -12 km/s and -50 km/s < U < -25 km/s), we get a 
mean disk metallicity of [Fe/H] - -0. 126 + 0.002 for the solar 
neighbourhood (the error representing the standard error of the 
mean). 

If we rather concentrate on the Hyades box, we get 
[Fe/H] = -0.061 +0.0 12 for Group 1, and [Fe/H] = -0.059 + 
0.013 for Group 2, meaning that the stream is over metal-rich 
at all masses. We can actually calculate the ratio of stars from 
the stream to stars from the field disk in the Hyades box from 
Figs. 3 and 4. We get A^h ya /Afield (box) = 3.65 for Group 1, and 
2.88 for Group 2, meaning that the stream has a roughly con- 
stant mean metallicity [Fe/H] = +0.02, even slightly rising 
for low-mass stars. This is inconsistent with the hypothesis of 
a field population from the solar neighbourhood trapped in the 
primordial cluster and subsequently evaporated. 

On the other hand, this constancy of the metallicity is con- 
sistent with the assumption that evaporated and field-like stars 
in the stream have roughly the same metallicity. This could im- 
ply that they both come from the inner regions of the Galaxy af- 
ter having crossed the corotation of the spiral pattern (the clus- 
ter could have been shifted in radius while remaining bound 
since the effect of a spiral wave on stars depends on the stars' 
phase with respect to the spiral, and the phase does not vary 
much across the cluster). The metallicity excess of 0.15 dex is 



indeed not consistent with the present-day orbit of the Hyades 
stream: this present-day orbit is centered on a guiding radius 
R g ~ 7.5 kpc in galactocentric coordinates (8 kpc being the 
galactocentric radius of the Sun), and would imply an implau- 
sible metallicity gradient of -0.3 dex/kpc in the disk. Assuming 
a very steep but realistic (see e.g. Daflon & Cunha 2004) galac- 
tic metallicity gradient of -0.07 dex/kpc, the metallicity excess 
is compatible with a galactocentric origin near R — 6 kpc. This 
would be consistent with the order of magnitude of stellar wan- 
dering caused by spiral perturbations (2-3 kpc, see Sellwood & 
Binney 2002; Lepine et al. 2003). A flatter metallicity gradient 
would imply an even more internal origin for the stream. 

5. Discussion 

Using the Geneva-Copenhagen catalogue of F and G dwarfs 
(Nordstrom et al. 2004), we have analyzed the mass function 
of the Hyades stream, an overdensity of stars in the f/V-plane 
kinematically associated with the Hyades cluster. 

In order to be compatible at a one-sigma level with the 
predicted mass function for stars evaporated from the primor- 
dial Hyades cluster, one needs a rather extreme scenario, with 
a very significant preferential evaporation from the cluster of 
1 M stars over 1.5 M stars, and with a very steep IMF, not 
flatter than a power-law exponent of -3.0 (see Sects. 4.1 and 
4.2). The initial total mass and virial radius of the cluster should 
also be very high (M > 1400 M Q , r > 7.5 pc). 

On the other hand, the observed mass function (see Sect. 
4.2, Eq. |9j is perfectly compatible with the hypothesis that 
the stream is uniquely composed of field-like stars obeying the 
PDMF. However, because the Hyades cluster is known to share 
the kinematics of the stream and to evaporate over time, the 
most likely scenario is that the stream is indeed mainly com- 
posed of field-like stars (about 85% of the stream for low-mass 
stars), but also partly of coeval stars evaporated from the pri- 
mordial Hyades cluster (about 15% of the stream for low-mass 
stars). 

The peculiar metallicity of the stream at all masses (see 
Sect. 4.3) is not compatible with a scenario where the field-like 
stars would have been trapped in the primordial Hyades clus- 
ter during its formation process, and subsequently evaporated. 
Those stars were thus most probably trapped at resonance by 
a spiral perturbation. Indeed, a series of strong transient spi- 
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rals with their mean corotation at the solar galactocentric ra- 
dius are known to produce small-scale structure in the local 
velocity distribution (De Simone et al. 2004). Stars on horse- 
shoe orbits that cross the corotation could wander over 2-3 kpc 
in much less than 1 Gyr (Sellwood & Binney 2002; Lepine et 
al. 2003), while staying on quasi-circular orbits, not elongated 
enough to betray their place of birth. This is consistent with a 
galactocentric origin near R = 6 kpc for the Hyades stream, 
that could explain its metallicity excess at low masses. On the 
other hand, the Hyades stream could also correspond to nearly 
closed orbits trapped at the 4 : 1 inner Lindblad resonance of 
a two-armed spiral density wave (Quillen & Minchev 2005). 
In any case, the prominence of the Hyades dynamical stream 
raises the question of the amplitude of the spiral perturbation 
needed to produce such a stream. The most detailed models of 
gas flows in the Galaxy (Bissantz et al. 2003) indicate that the 
amplitude of the spiral structure in the mass density is larger by 
a factor 1 .5 than its amplitude in the near-infrared luminosity 
density, while such a large amplitude is also needed to produce 
the large non-axisymmetric motion of the star-forming region 
W30H (Xu et al. 2006). This requires the baryonic disk to be 
massive even near the Sun, and is very constraining for the dark 
matter distribution in the Galaxy (see e.g. Famaey & Binney 
2005). On the other hand, the probable resonant origin of the 
Hyades stream cautions against naive backward integration of 
individual orbits in an axisymmetric potential when describing 
the past evolution of such streams, and it bares the importance 
of evaluating the impact of radial migrations on the chemical 
evolution of the Galaxy (see also Haywood 2006). 
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